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Reflection of electromagnetic waves from a stratified lossy inhomogeneous plasma is discussed. 
The profile for the conductivity is idealized as an exponential curve with a superimposed Gaussian- 
shaped perturbation. By examining the change in the calculated reflection coefficient for various 
locations of the perturbation, some insight into the reflection process is gained. In particular, it is 
shown that reflection does not take place at a single level but, instead, a wide range of levels is im- 
portant. The parameters of the problem are chosen to be representative of the D layer of the ionos- 
phere and the wavelengths cover the range 10 km to 30 km. 



1. Introduction 

It is well known that the reflection coefficient of a 
stratified medium is determined explicitly by the profile 
of the complex dielectric constant. For example, if 
the profiles of electron density and collision frequency 
in the lower ionosphere are specified, it is possible to 
calculate the reflection coefficient for an obliquely 
incident electromagnetic wave. Although there is 
now a vast literature devoted to this general subject, 
some questions concerning the nature of the reflection 
process remain unanswered. For example, the deter- 
mination of the effective reflection level in the ionos- 
phere for VLF radio waves is not a clear-cut task. In 
fact, even the concept that waves are reflected at some 
specific level is open to question. 

It is the purpose of the present paper to investigate 
which levels in the lower ionosphere are most influen- 
tial in determining the reflection coefficient. A num- 
ber of simplifications are made in order that extraneous 
considerations are removed. For example, in the 
lowest ionosphere in the daytime it is assumed that the 
following inequalities hold: 

(a) Electron collision frequency v e > electron 
gyrofrequency ft e - 

(b) Ion collision frequency v\ > angular wave fre- 
quency 0). 

(c) Density of neutral particles N > both density 
of electrons N e and density of ions Ni. 

Under these conditions, it is then a straightforward 
matter [Wait, 1962] to show that the ionized medium 
may be characterized by a scalar complex dielectric 



constant e(z) which is a function of height z of the form 
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In the above, 



€o is the dielectric constant of free space, 
e is the charge of the electron and the ion, and 
m e and m\ are the masses of the electrons and the 
ions, respectively. 

The only difference between this formulation and 
a previous one [Wait and Walters, 1963], is that here 
the ions are considered explicitly. The form of the 
dielectric constant €(2) is just the same. Thus, the 
numerical method for calculating the reflection coeffi- 
cient is identical and need not be described here. 

2. The Idealized Profile 

To illustrate the salient features, it is assumed that 
the conductivity parameter l/L(z) of the ZMayer may 
be written in the following analytical form: 
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FIGURE 1. The idealized profile of the inhomogeneous lossy plasma. 

The ordinate is proportional to the conductivity of the medium as a function of height 
above and below the reference level at 2 = 0. The perturbation is shown at a number of 
levels from z = F = — 8 km to +4lkm. 
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The resulting exponential profile obtained by setting 
A = may be described as the ''undisturbed profile." 
Typically, for daytime conditions, /3 would be about 
0.3 km -1 when z is measured in km above some refer- 
ence level at z = 0. The second term in square brack- 
ets is then regarded as a perturbation which is centered 
at 2 = F, while its vertical extent is determined by the 
parameter D. In an abnormally ionized D region, such 
as discussed by Crain and Booker [1964], the ions can 
contribute significantly to the perturbation. The 
function l/L(z), as defined above, is plotted in figure 
1 as a function of height z where Lq is set equal to 1/2. 
The smooth curve in this figure corresponds to the 
undisturbed exponential profile. The effect of 
locating the perturbation at various heights is illus- 
trated in figure 1 by choosing a sequence of F values 
from — 8 km to -f 4 km, while D is chosen to be 0.5 km 
and A is taken to be 0.5. 

It is stressed that the perturbation, as described 
above, is highly idealized. It is not intended that it 
should represent the Z)-layer under actual disturbed 
conditions. 
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Figures 2a and 2b. The (vertically polarized) reflection coeffi- 
cient as a function of F which is the vertical level of the perturba- 
tion. 



3. Description of Numerical Data 

It is now contended that considerable insight into 
the reflection process is obtained by examining the 
changes in the reflection coefficient as the perturbation 
is moved along the profile as indicated in figure 1. 
Results of this kind are shown in figures 2a and 2b 
for the amplitude and phase of the reflection R for 
a vertically polarized incident wave. The cosine C 
of the angle of incidence is denoted on the curves. 



Also, as indicated, the wavelength A is taken to be 
15 km corresponding to a frequency of 20 kc/s. The 
results for the undisturbed exponential profile are 
shown in figures 2a and 2b by dashed lines. Follow- 
ing previous conventions [Wait and Walters, 1963], 
the phase is referred to the level 2 = 0. 

It is evident from the curves in figure 2a that the 
perturbation only has an influence on \R\ if it is located 
below the level 2 = F = 0. For F in the interval from 
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about — 1 to — 6, the reflection coefficient is increased 
somewhat. Thus, the perturbation at these levels 
actually enhances the reflection process. Of course, 
for increasingly larger negative values of F, the 
opposite tendency becomes evident. Then the pertur- 
bation acts as an absorber, tending to diminish the 
reflection coefficient. 

It is suggested that the "levels of reflection" in 
the unperturbed profile are characterized by the 
interval in F where | R | for the perturbed profile ex- 
ceeds \R\ for the unperturbed profile by some detect- 
able increment. For example, when C = 0.2, the 
"reflection levels" extend from about F = —2 to 
— 10 km. In this range, \R\ for the perturbed profile 
exceeds \R\ for the undisturbed profile by at least 
5 percent. It is interesting to compare this deduction 
with a statement of Crain and Booker [1964] that the 
reflection takes place in the height range near the level 
given by 
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For the undisturbed exponential profile, this is equi- 
valent to solving the equation 

2 exp (fiF) = 2C 2 . 

Thus, taking (3 = 0.3 and C = 0.2 as in the above ex- 
ample, it is readily found that 

F=-ll km. 



waveguide propagation of VLF radio waves is the 
complex quantity a [e.g., Wait and Spies, 1964]. In 
terms of the reflection coefficient R, it is defined by 

R=— exp (aC), 

where C is the cosine of the angle of incidence. 
Writing a = cti + ia2, the real quantities ot\ and «2 
are shown in figures 3a and 3b for C = 0.l and 0.2, 
where we have used the perturbed exponential model 
described above. The magnitude of ot\ is approxi- 
mately proportional to the attenuation of the domi- 
nant mode in the earth-ionosphere waveguide. Thus, 
for F greater than zero, the attenuation is not influ- 
enced by the perturbation. However, for a range of 
F from about — 1 to — 10 or so, there is a detectable 
change in the attenuation rate. 

The variation of the reflection parameters a\ and 
«2 as a function of the angle of incidence is shown in 
figures 4a and 4b for a range of F values. It is im- 
portant to observe that the functional dependence on 
C is markedly changed when the perturbation is 
lowered from F = to F = — 10 km. The consequence 
is that the relative attenuation of the first and second 
order waveguide modes is appreciably modified by a 
low-lying perturbation of a smooth undisturbed profile. 
The general tendency is to increase the relative attenu- 
ation of the second order mode as the perturbation is 
lowered. The phase parameter cti in figure 4b indi- 
cates that the perturbation generally increases the 
phase velocity of the dominant waveguide modes. 



Crain and Booker [1964] do not explain how they arrive 
at their criterion for the reflection level but, pre- 
sumably, it is related to the phase integral approach 
of Eckersley [1932]. The latter is straightforward for 
lossless media at high frequency where such geo- 
metrical-optical notions are meaningful. For lower 
frequencies in relatively lossy media, the phase integral 
method loses its validity. This would appear to be the 
reason for the apparent discrepancy between the 
Crain and Booker "reflection level" at F=— 11 km 
and the results in figure 2a which would indicate that 
reflection takes place at a range of levels from about 
— 2 km to — 10 km. The difference between the phase 
integral approach and the full wave solution becomes 
even greater for angles nearer grazing incidence. 
In general, the reflection level predicted by the Crain 
and Booker criterion is at the bottom or somewhat 
below the range of levels in the lower ionosphere where 
VLF radio waves are reflected. However, this view 
does not modify the conclusions regarding the import- 
ance of ions which was the principal topic in the paper 
by Crain and Booker [1964]. Furthermore, their 
"reflection level" is probably quite representative at 
the height region where the bulk of the energy is act- 
ually returned. 



4. Connection With Mode Theory 

An important parameter occurring in the theory of 
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FIGURES 3a and 3b. The re/lection parameters ot\ and «2 which are 
defined by R= — exp (aC) where a — ai -fi ol<l. 
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FIGURES 4a and 4b. The reflection parameters as a function of the 
cosine of the angle of incidence. 



For example, at a frequency of 20 kc/s, the phase 
velocity (relative to c) of mode n=l changes from 
0.9986 to 0.9989 as the location of the perturbation 
moves from F = to —10 km. In this case, it is 
assumed that the C = 0.15 and the reference level 
(z = F= z 0) in the ionosphere is located at a height 
h = 70 km above a perfectly conducting earth. Other 
quantitative estimates of this kind may be obtained 
by using the waveguide parameter curves issued in a 
recent NBS publication [Wait and Spies, 1964]. 



5. Wavelength Dependence 

The results discussed above all pertain to a wave- 



length of 15 km. To illustrate the wavelength depend- 
ence, the amplitude and phase of the reflection coeffi- 
cient are plotted as a function of C in figures 5a and 
5b, respectively, for wavelengths extending from k= 10 
km to 30 km. The function L is defined by L = 30/X 
where \ is the wavelength in km. It is interesting to 
observe that the curve in figure 5a for k = 10 km has 
reflection levels extending from F = —2 km to well 
below F = — 10 km. The Crain-Booker [1964] reflec- 
tion level would be found by solving 

3exp08F) = 2C 2 = O.O2, 

which yields F = — 16.5 km which is not far below the 
physically important height range where the reflection 
actually takes place. 
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FIGURES 5a and 5b. The wavelength dependence of the reflection 
coefficient. 
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6. Conclusions 

The graphical results given in this paper would 
indicate that reflection does not take place at a dis- 
crete level in an inhomogeneous lossy medium with a 
smooth profile of its dielectric constant. In fact, 
the wave is reflected at a range of levels in the medium. 
The upper limit is determined where the incident wave 
has been sufficiently attenuated by a result of conduc- 
tion losses. This is in contrast to the mechanism 
which produces reflection in a lossless inhomogeneous 
medium. In that case, there is a fairly well deter- 
mined reflection level which is located where the 
vertical propagation constant vanishes. 

It is important that the reflection process be under- 
stood if experimental data on VLF fields are to be 
interpreted in relation to the properties of the lower 
ionosphere. All too often, variations of the amplitude 
and the phase of the VLF field are ascribed to changes 
in the parameters of some simplified model which 
bear little resemblance to the actual physical changes 
which have taken place. 
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